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Abstract 

Analytical, free of time consuming Monte Carlo simulations, framework for credit portfolio sys- 
tematic risk metrics calculations is presented. Techniques are described that allow calculation of 
("^ ' portfolio-level systematic risk measures (standard deviation, VaR and Expected Shortfall) as well as 

CN| , allocation of risk down to individual transactions. The underlying model is the industry standard 

multi-factor Merton-type model with arbitrary valuation function at horizon (in contrast to the sim- 
plistic default-only case). High accuracy of the proposed analytical technique is demonstrated by 
benchmarking against Monte Carlo simulations. 



1 Introduction 



There exists an increasing demand for fast and consistent economic capital calculation and allocation 
Cr , • techniques. Portfolio-wide calculations of economic capital are just a first step in the modern process of 

credit portfolio management. Financial institutions are more and more involved in stress testing, sensi- 
tivity and scenario analysis. For these purposes the portfolio-level risk measures need to be recalculated 

over and over again. Using industry standard Monte Carlo simulations for the portfolio-level risk quantifi- 

^^- cation requires considerable amount of time and computer power. For the purposes of risk concentration 

identification, risk-adjusted pricing and portfolio optimization the portfolio-wide risk (economic capital) 
needs to be allocated down to individual transactions. The latter task is even more challenging from 
^ I both methodological and computational points of view. Statistical noise, being inherent part of Monte 

Carlo simulations, leads to unstable estimations of the allocated risk (especially in case of VaR-bascd cap- 
ital allocation). Reliable estimations of capital charges based on simulations require significantly more 



iy-\ ■ computer time/power compared to the portfolio-wide calculations. 



Although several techniques have been developed t o improve the pe r forma nce of the simulations- 
based ap proach, e.g. importance sampling (see, e.g., iKalkbrener et al.L 120041) and kernel estimators 



f^ I (see, e.g., iTaschel 120091 ). simulation-based estimation of risk contributions on transaction level is still a 



demanding computational problem. In practice, when applied to multi-factor models, efficiency of the 
abovementioned techniques may be limited since it depends on a quality of analytical approximation used 
- to determine a sampling region. Yet another drawback of the simulation-based approach is its inability 

/\ ' to efficiently risk-assess new deals in a context of the portfolio. 

S . An alternative to the simulation-based approach would be some kind of analytical technique. Although 

Merton-type models are not analytically tractable in general case, some progress has been made to develop 
an approximate solution. Th e most s uccess ful attempts to tackle the problem are Asymptotic Sing l e Risk 
Fac tor (ASR F) fram ework (JGordvl . [20031) . granularity adjustment (GA) by iMartin and Wildj (|2002[) 



and iPvkhtinl 's ( 2004 ) m,ulti-factor adjustment. This article aims to complement the existing analytical 
techniques. The ambition is to fill the existing gap betwe en theoretical results and practice by considering 
a fully- featured PortfolioManager-type (JKealhofen . 120011) credit portfolio model. The proposed framework 



allows to calculate most commonly used risk measures (variance, value-at-risk and expected shortfall) 
on both portfolio and transaction levels. The material presented here lacks mathematical rigor. Instead, 
results of numerical tests are presented to demonstrate the performance and prove the validity of the 
proposed techniques. 



'Independent consultant. E-mail: |mvoropaev@mail.ru| 



This article is organized as follows. First, in Section[2l a short description of the multi-factor Merton- 
type model is given, followed by a review of the progress made so far on the model's analytical tractability. 
VaR expansion technique, used as a starting point for the approach presented here, is presented in Section 
[31 Main results are presented in Section 21 where conditional expectation scries expansion is derived 
(Sections 14. II and 14. 2|) and applied to systematic risk constituents. It is demonstrated, how the proposed 
expansion technique can be utilized to compute systematic components of various portfolio-wide risk 
measures and corresponding risk contributions. Finally, in Section [5l the conditional expectation scries 
expansion technique is extended to cover idiosyncratic risk components. Monte Carlo simulations are 
used to substantiate the validity of the proposed analytical approach (Sections 14.51 and 15. 4p . 

2 Structural credit portfolio models 

Merton-type credit portfolio models are most widely accepted ones for the purposes of credit portfolio 
risk metrics calculations. In these models the portfolio consists of risky instruments {vi} with the value 
Vi of each instrument at horizon (usually set to one year) being a function of normally distributed random 
variable e^ (normalized asset return). Correlations between these variables {ci} are modeled through a 
set of Nf normally distributed independent variables {rjk} referred to as common factors. Each variable 
Si is split in a sum of instrument specific (idiosyncratic) part, which depends on a Gaussian variable ^i, 
and systematic part, which depends on the common factors, as follows 



v^{e,) = z;.(p^Efc(/3.)fc^/fc + V^ - Pf 6)- (2-1) 

The independently distributed random variables {{Ci}? {'?fc}}u a-re assumed to have zero mean and unit 
variance. Instrument specific constants |pi| < 1 and {{Pi)k} determine dependency of e^ on the common 
factors (related to geographic regions and industry types). The so-called factor loadings {(A)fc} are 
subject to normalization condition 

EJft)?- = 1- (2.2) 

Uncertainty in the value of the portfolio V = ^^ Vi is quantified by means of various risk measures, 
most popular of which are VaR(Value-at-Risk) , ES(Expected Shortfall) and standard deviatioro. 

Once the portfolio-level risk measure is known, the question arises how to distribu te (allo c ate) t his 
risk consistently among the constituents. The Eulcr allocation technique (see, e.g., iTaschd . 12008 ) is 



the commonly adopted solution. According to the Euler allocation principle, individual assets Vi of the 
portfolio are assigned fractions (risk contributions) 9i of the portfolio-level risk 6 according to 



96 

dwi 



e = ^0„ (2.3) 



where Wi is a weight of ith facility in the portfolio. In what follows the weights {wi} will be implied but 
not written explicitly. 

No closed-form solution exists for either portfolio-level or facility-level risk measures in the general case. 
Several important steps have been made towards approximate analytical solution of the problem. First, 
the case of one common factor and infinitely large and fine-grained portfolio was solved by Asymptotic 



Single Risk Factor framewo rk (iGordvl. l2003f). Next, i diosyncratic component of risk has been ad d ressed 



by granularity adjustment ( Martin and Wildel. [2002 ) . Finally, the results of iMartin and Wildd ( 20021) 



were applied to a multi- factor case bv IPvkhtinI (2004) 



Unfortunately, no significant progress has been made ever since towards better a nalytical approxima- 
tion; h owe ver, some generalizations of previous results have been recently reported bv lGordv and Marrond 
(|2010l ) and lOagliardini and Gourierouxl (120101) and attempts have bee n made to find a more simple solu- 



tion to the multi-factor case bv lDuellmann and MasscheleinI ( 20061 ) and lCesoedes et al.l ( 200(31) . Moreover 



practitioners considering applying Fykhtin's approach to realistic credit portfolio models face two major 



^Assuming {^i} to be independently distributed is equivalent to an assumption that eaeh borrower in the portfolio is 
represented by one faeility. This assumption is made to simplify notations and does not undermine the validity of the 
results. 

^ See iHuIl ||2007| ) for a detailed diseussion of the various risk measures. 



difficulties. First, Pykhtin's model was formulated for a default-only case and it is not at all obvious 
how to (efficiently) extend it to a more general (and realistic) case of value-based valuation at horizon. 
Second, calculation of the multi-factor adjustment are of quadratic in portfolio size complexity, making 
application of the model to large portfolios barely possible. On top of that, no solution to the problem 
of risk allocation within Pykhtin's model has ever been reported. 

In the following sections a new approach is presented. Although based on the same principles, the 
approach will address the above mentioned difficulties of Pykhtin's model. Moreover, higher order (i.e. 
third order vs. original second order) multi-factor adjustments will be considered. Many of the restrictive 
simplifications (e.g. default-only mode, simplistic correlation structure, homogeneous portfolio, etc.) 
commonly used in previous publications will be loosened, allowing application of the proposed analytical 
framework to a wide class of realistiqfl structural credit portfolio models. 

3 VaR and ES adjustments 



Building on the work of ICourieroux et al.l (|200CI| ). iMartin and Wildd (J2002I ) derived the second order 



correction to VaR and used the results in the context of credit portfolio to calculate an adjustment for 
undiversified idiosyncratic risk {granularity adjustment). Somewhat simpler derivation is presented here, 
outcome of which is a higher precision correction to VaR and is more suitable for the techniques presented 
in this article. 

Consider random variable x with continuous probability distribution function (p.d.f.) f{x). Let qa 
be the a-quantile of this distribution. Consider another random variable Sx with g{Sx\x) being its p.d.f. 
conditional on the value of the first variable x. Let us find the a-quantilc (7* of the p.d.f. f*{x + 6x) of 
the sum of the above two variables. The /* can be written as 



/*(-^) = / fix - 6x)g{dx\x - Sx)d{Sx) (3.1) 

Expanding the right hand side of this expression in Taylor series of {x — Sx) around x, one can obtain 

fix) = fix) + Y, -^^^f^''^^'"(''^^^ ^"(^•) = / {Sxr9{Sx\x)d{Sx), (3.2) 

n— 1 

where ^n{x) are moments of 5x distribution conditional on x. 

Once the relationship (j3.2p between probability distribution functions has been established, the rela- 
tionship between quantiles can be derived by substituting p.2p into the following definition of a-quantile 

"' f{x)dx^ T" r{x)dx (3.3) 

— 00 J —OQ 

The result is 

ril '^ ('— Ti"~i d"~i 



gn 



2^ — 9a 



Suppose Sx is a small correction to x. One way to quantify this smallness is to assume that /i„ ~ ^", 
where S is some small number. One can solve the equation p.4p order by order in S by expanding its 
both sides in powers of (q* — Qa) around qa- 

Only distributions satisfying fJ-i{x) = will be considered in this article. In this case the {/i„(x)} 
become conditional central moments and p.4[) has a particularly simple third order solution 



d 



2/(a;)da;''" "^"' " x=q^ 6f{x)dx^ 



2 



9a-9a«-77777T3z[/(2;)M2(x)] +F77ZTT::7[/(^V3(a:)] (3.5) 



X — Qa 



Let us look at the result p.5[) from credit portfolio perspective. Let a; be a single factor approximation 
of the portfolio value, x = V{t]). Let the factor rj be normally distributed with the p.d.f. n{ri) = 



'realistic" here means "used in practice" 



e ^ /^/\/27r. The a-quantile qa is related to the portfolio's VaR and portfolio's expected value E{V) aQ 

VaR = E{V) - Qa (3.6) 

Using n'M = —rin{r]), f{V)(W = n{r])dri and (|3.5p . the second and third order VaR adjustments can 
be written aa 



AVaR2(a) = 



AVaR3(a) 



1 d fn{ri)fi2{ri) 



2 n(ry) dr/ \ V'(ri) 
1 / , / V" 



7]=^^^{a) 



i;=$^i(a) 



1 



1 d fn{T])n3{T]) 



6n{ri)dri \V'{r]) drj \ V'{ri) 



6\V 



V 



2 



1 f,J_,;f2, + 3X:),„3f(-,-l) + 3X + ?in-'"^"' 



v 



[V] 



(3.7) 



(3.8) 



?7=<1>^1(q;) 



where ^~^{a) is the inverse of the normal cumulative p.d.f. 

Using the VaR adjustments (|3.7p and p.7p . one can easily calculate similar adjustments to expected 
shortfall. Noticing that 

ES(a) = - / YaR{ri)n{rj)dri, (3.9) 



(3.10) 
(3.11) 



the second and third order expected shortfall contributions can be written as 

1 n 



AES2(a) 
AES3(a) 

4 Systematic risk 



H2 

2aV' r,=*-i(Q) 

In/, / V" 



?)=<E>~i(a) 



Let us start by assuming that the portfolio dynamics is mainly governed by the systematic risk compo- 
nents, i.e. the common factors {rjk} give the main contribution to the portfolio risk measures, while the 
idiosyncratic factors {^i} give rise to less significant corrections. In this section it will be demonstrated, 
how to isolate the systematic risk by integrating out the idiosyncratic components. It will also be shown, 
how the resulting asymptotic multi-factor framework can be utilized for the risk metrics calculations. 

The same framework will be further extended in the next section to cover the idiosyncratic risk 
components. 

4.1 Series expansion for conditional expectation: single factor 

In order to focus on the systematic part of portfolio dynamics, let us integrate out (average over) the 
idiosyncratic component ^^ in (j2.1|) . Let us assume there is just one common factor and extend the results 
to a multi-factor case later. 

Average value of a facility ui conditional on the systematic factor rj is 



^>^{v)= hiiPiV+ V^ - pH)- 



-e/2 



/27r 



-de, 



(4.1) 



*VaR defined this way is simply an economic capital of the portfoHo. 

^The signs of both VaR and ES adjustments seem different from those that can be found in the Hterature. This apparent 
contradiction is explained by the fact that the analysis here is based on the value of the portfolio V, rather than its losses. 



which after changing the integration variable to e = piri + -^/l — pf ^ becomes 

1 /2p,ery-pf(e2+,y2)\ ^-e^/2 



Viiv) = v{e) 



V^ 



: exp 



2(1 -Pf) 



/27r 



de. 



(4.2) 



The above e xpression can be further developed by applying Mehler's formula (for the proof see, e.g., 
Foatalll978[) : 



^He„(e)He„(77)^ 



n=0 



V^ 



zcxp 



2(1 -p2) 



(4.3) 



where He„(77) = ( — l)"e'' ^■^(d/dri)'^e ^ 1"^ are Hermite po lynomials (for definition and properties of Her- 
mite polinomials see, e.g.. lAbramowitz and Stegunl . ll972i) . The result is 






,(") 



Vj(e)He„(e) — ^^de 



/27r 



(4.4) 



A few remarks are needed regarding the result (j4.4p . First, expansion exists as long as all the coef- 
ficients u^") are finite. This is the case, for example, for any piece- wise continuous function Wi(e) whose 
absolute value at infinity (e — >■ ±oo) does not increase faster than some power of e. Any reasonable value 
function of a financial instrument does satisfy this constrain. 

Next, in case p = 1, the classical Hermite series expansion is recovered. The series converges to 
the value of the function everywhere except for discontinuity points where the series converges to the 
average of the function around the discontinuity point. The Hermite series expansion is known for its 
slow convergence especially for large values of the argument 77. 

Finally, as a consequence of |p| < 1 in (|4.4[) . the conditional expectation series converge significantly 
better. For the same reason, i.e. \p\ < 1, the conditional expectation function uiir]) is not only continuous, 
but differentiable infinite number of times. 

Before generalizing the result (|4.4p to a multi-factor case, let us explore the benefits of the expansion 
(|4.4p in the context of credit portfolio. Advantages of the proposed approach can be seen even in a simple 
case of a single factor model. 

The asymptotic single risk factor 77 value of the portfolio ^1/(77) — ^^ W(^/) can be easily derived from 
(|n)) and is 

Vifiv) = ^F(")He„(r;), F^ = ^ ^t;^) (4.5) 

n i 

Once the coefficients F'-"-' are calculated, one can immediately write both VaR and ES of the portfolio 
for any confidence level a as 



VaR(a) = -^F(")He„(77) 



r)=<E>^i(Q) 



ES(a) 



.-vV2 
V2^ 



5^y(")Hc„_i(r7) 



i]=i~'^{a) 



(4.6) 



Using 



ri>0 / V / V "" „^g 

and (|4.5p . trivial calculations lead to the following VaR and ES -based risk contributions 



VaR,^ = -^^T;(")He„(r7) 



n>0 



JJ = $^l(C() 



Es; 



-r,V2 



(") = ^^EV"^Hc„_i(77) 



n>0 



?)=<E>^i(a) 



(4.7) 



4.2 Series expansion for conditional expectation: multiple factors 

In a multi-factor case, the conditional expectation (|4.4p can be written as 



W('7/.) = E^^'"^Hc„(E,(/30fc'7fc) 



(4.8) 



This expression, however, does not allow to write the portfolio value V in a form similar to ()4.5|) . To 
accomplish this, let us introduce multivariate Hermite polynomials 

n 

The multi-factor expansion then becomes 

W(%) = E 4^^^"^ (^^)*^^(Sr"He^^2 . . .(,^^,) (4.10) 

n 

and the conditional expectation of the portfolio can be written as 



n..)-E^^E <L_,^«" "•■(^)' ^^-E^-f^(/30..(A)..... (4.11) 



Using orthogonality properties of multivariate Hermite polynomials 



J He^^-- fa)He^'-- (7?fc) ^ ^^ J^;' J^ drjk = n! 5„„,4,/,<5,,z, . . . , (4.12) 

one can calculate the variance (Jy of the portfolio 

a^ = E(y^)-(E(y))^ = E-! E [<LX (4-^3) 

n>0 /£i,fe2,... 

Standard deviation ay based risk contributions can be calculated using the (|2.3|) and (|4.1ip . The result 

is 

^.— -^E/^r^i"^ E (A:).. (A)..... <!,... (4-14) 

av — — 

n>0 kiM,--- 



Recently, it was shown bv lVoropaevI (|2009l ) that applying (j4.13|) and (|4.14p results in calculations which 
are of linear complexity in portfolio size. The amount of common factors Nf, however, is the bottleneck 
of the calculations. Indeed, nth term in the above expressions contains N^ elements, making calculations 
of higher order terms impractical. Fortunately, only a few first terms lead to an acc urate resul t s. Fo r 
details and discussion of the convergence properties of (|4.13p the reader is referred to IVoropaevI ( 20091 ) . 



where the problem of standard deviation and standard deviation base risk allocation has been solved in 
more general case using techniques similar to those described here. From now on we will focus on the 
tail risk measures, VaR and ES. 

4.3 Conditional expectation in the tail 

Let us assume that the portfolio value distribution in the multi-factor case can be approximated by some 
single-factor value distribution, i.e. let us write the value of the portfolio as 

V = V^f{Y) + Vmf, E(K„/|Fi/)=0, (4.15) 

where Vif is a single-factor approximation and Vmf is a multi-factor correction with zero expectation 
conditional on Vif. The single systematic risk factor F is a linear combination of the common factors 
{r^fc}. The choice of the principal risk factor Y is somehow arbitrary; however, one would aim to choose 
Y such that Vif is as good approximation to V as possible and Vmf is as small correction as possible. A 
solution to this optimization problem (which needs to be well formulated first) may be a matter of future 



research. Fortunately, as we will see later, even in case of sub-optimal choice of Y one can achieve very 
good numerical results. 

The (sub-optimal) choice of Y used here is based on the following rationale. Notice that nth term in 
the conditional expectation expansion (j4.1ip is (roughly speaking) proportional to p" where pp is some 
characteristic correlation. Assuming pp is small, one can conclude that the lower order terms in (j4.1ip 
give the main contribution to the portfolio dynamics. Assuming further that the n = 1 term is the most 
important one, one would naturally choose Y to point in the direction defined by 

VW = ivl'\v^'\...,V^';) (4.16) 

This particular choice of Y is not only natural and convenient within the proposed framework, but also 
will be substantiated by numerical tests in Section H?gFl . 

One last preparation is needed before splitting the portfolio value according to (|4.15p . Once the 
principal risk factor Y is known, let us transform the initial orthonormal set of common factors {rjk} 
by some orthogonal transformation in such way, that one of the transformed factors coincides with Y. 
This can be achieved by Gram-Schmidt process starting with Y. From now on we will assume that the 
transformation took place and that {rjk} is a set of the transformed common factors. The rji factor is 
assumed to be the principal risk factor. 

To split the portfolio value (|4.1ip according to (|4.15p . let us make use of the following identity, which 
can be derived using the definition of the multivariate Hermite polynomials (|4.9p and the fact that V^ "j! 
are symmetric in ki,k2, ■ ■ ■, 

n 

y^:l.kM'-'"M=T.^n. Ki.. (';)He„_,(^i)He,(,yD, (4.17) 

1=0 

where (") are binomial coefficients and rj'^ is a set of all common factors but rji. Using the above 
expression, the portfolio value (|4.11[) can be written as 

Y^ Y^ Y^ /^""^w. ^.,M/(™) VfMk2 




n ki,k-2,...m>n 



V{v.)=Y: E E „ H^™-('nKr..fc,fc,...He^"---(^D. (4.18) 



Finally, separating the n = term and introducing conditional coefficients V^"J {vi)i the portfolio 

value can be put into the form 

ViVk)^V,firji) + V^M\m) (4.19) 



^i/(^/i)=E^i/Hc„(^i), K™/(r?:-hi)-E E ^rnf\^k.JVi)iiot^''-iril) (4.20) 

n n>0 ki,k2,... 



m>n 



The multi-factor correction Vmf in the above has zero expectation conditional on 771 due to the orthogo- 
nality properties (|4.12[) . For a given confidence level a, the above expressions represent series expansion 
of the conditional (on 771 ~ $~^(a)) tail expectation. 



^Pvkhtir] I2OOJ) suggests different choices for Y . These choices, however, are intuition-based and are not theoretically 
substantiated. In the author's experience, the choice of Y presented here leads to better results when applied to realistic 
portfolios. 



4.4 Systematic tail risk and its allocation 

The series expansion of the conditional tail expectation (|4.19p - (|4.2ip together with the single- factor case 
results (|4.6p allow us to apply the results of Section [3] to VaR and ES calculations. 

Since the single- factor VaR and ES have been calculated before, i.e. (|4.6p . let us start with the second 
order contributions p.7p and p.lOp . Using the notations introduced in the previous section, the second 
order VaR and ES adjustments are 



AVaR2(a) 

AES2(a) 



1 



2^i'/(^i) 
1 n{rii) 



^'2ivl)-^^2{vl)U+^^^ 



T;i='l'^i(a) 



M2(?7l 



7,1 =*- 1(a) 



The Vif derivatives can be calculated using (|4.20p and are 

^Ifim) = E vlfnRe^_,{in), Vl'fim) = Y. Vlfn{n - l)He„_2(77i) 

n>0 n>l 

The conditional second central moment (variance) ^2('/i) is 

n>0 fci,fc2,... 

and its derivative fJ-2ivi) can be calculated as 



(4.22) 



(4.23) 



(4.24) 



M^(m) = 2E-! E eA,.,...(^i) 

n>0 ki,k2,... 



V 



(n) 



™/fclfc2... 



(^71) 



(4.25) 



where 



v: 



(») 



™/ fcifca 



im 



E 



(m) 



(m - n)Re,n-n-iim)Vii...k,k; 



(4.26) 



The above solves the problem of second order VaR and ES adjustments on portfolio level. The 
corresponding risk contributions can be calculated by applying (|2.3p to (|4.22p . This exercise is left for 
the reader who may find useful the following examples 



d 



^'7j;7^i/('?i) - Err"^ (^^)>Hc„-i(ryi) 



' dw; 
d 



n>0 



(4.27) 



— A.2(r;i) =2Epr E Ki"A,,,...('?i)E LjHc— (m>^^(ft)r"(/30..(A).: 



n>0 ki,k2,... 



m > n 



Calculations of the third order VaR and ES adjustments, (|3.8p and (j3.1ip . and corresponding risk 
contributions can be done in the same fashion. The difficulty one will face in this case is calculation of 
fi^^rji). To calculate the third central moment the following integral has to be evaluated 



M3(??l 



ynifiVklm) 



(2^)(^/-i)/: 



-dr/fc 



(4.28) 



Unlike the case of /^2(??i), orthogonality conditions (|4.12p alone are not sufficient to calculate the integral. 
One is facing the problem of calculating exponentially weighted average of three Hermite polynomials. 
To solve this problem, let us start with the following identity (which follows from a more general result 

fc!He„+™_2fc(a;) 



io solve this pr 
of [PrakQ (i2009h 



He„(a;)Hem(.T) = E 



k)[ k 



(4.29) 



The integral then can be solved as follows 



e ^' 1"^ n\vti\k\ 



da;He„(.T) Hem(.T) Hefefx) — ^^ = -. — —, r — 7^- r — -. — ; r^, (4.30) 



2 y V 2 y ■ V 2 



provided m + n + k is even and each oi m,n,k does not exceed the sum and is not less than the absolute 
value of the other two. Otherwise, the integral is zero. 

It is not clear how to write multivariate version of the above identities. However, using the above 
identities together with the definition of multivariate Hermitc polynomials (j4.9p . one can solve for any 
given set of n, m, k. For example, 

He^ ivm4''' ivDM'''"' iv*k) l2n)iN^-i)/2 d^fe = QSk,kJk..k.Jk,k, (4.31) 

First few terms of the third central moment /i3 (rji ) are 

f^sM = 2EC(^i)<\.(^i)^™/.,..('?i) (4-32) 

ki,k2 

+ 8 E ^k.k.(:^^)^k.k.^^)^k.kM) + --- 

ki,k-2,k3 

The results presented in this section allow to calculate portfolio-level and facility-level systematic 
components of VaR and ES. It is easy to see that the necessary amount of calculations is linear in a 
number of facilities of the portfolio. Moreover, the calculations can easily be parallelized on a multi- 
processor machines. 

4.5 Numerical results 

To prove the validity and demonstrate the accuracy of the proposed analytical framework, let us compare 
results of the analytical approximation with those of unbiased Monte Carlo simulation. The focus here 
will be on VaR and VaR-based risk contributions. Only the systematic risk component is considered here. 
Numerical tests covering idiosyncratic risk as well will be presented in Section 15.41 

Since we are interested in the systematic components of portfolio risk, the Monte Carlo routine used 
here was developed to cover systematic, but not idiosyncratic risk components. This is achieved as 
follows. For each Monte Carlo scenario a set of systematic factors is generated. Instead of generating 
borrower-specific factors, however, expected (given systematic factors) values are assigned per facility. 



The particular set of common factors used in the tests is similar to the one described in iKealhofer 



(J2001I ). The total of Nf ~ 120 factors cover 61 industry and 45 regional sectors. Two portfolios were 
constructed, diversified and concentrated. Both portfolios contain identical loans maturing at horizon. 
Each loan's correlation with the systematic factors pi is 0.6 and probability of default (PD) equal 1%. 
The corresponding value function z'i(e) is 

, , f 1 if e> $-1(0.01) ,, ^^, 

<'^-[o if. < $-(0.01) (4-33) 

The diversified portfolio contains 45 x 61 = 2745 loans, each loan representing a different region/industry. 
The concentrated portfolio contains 400 loans randomly assigned to different region/industry and 100 
loans representing a single region/industry pair. These 100 loans create region/industry concentration in 
the portfolio. 

Monte Carlo estimates of portfolio VaR and VaR contributions per facility were based on 10^ scenarios. 
Confidence interval was set to 99.9%. Estimates of VaR contributions were calculated based on 50,000 
scenarios around the 99.9% point (i.e. average VaR contributions for 99.875% - 99.925% interval were 
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Figure 1: Relative differences between Monte Carlo and analytical esti- 
mates of the systematic VaR-based risk contributions. 



calculated). Plain vanilla Monte Carlo simulations were used to exclude any bias and limit possibilities 
of implementation errors. 

Several analytical estimates were calculated. First, single factor approximation (If) was calculated 
based on (|i?TO)) - (|T^ and (|i^ - (|i?7)) . Next, second order (multi-factor) VaR adjustment (|X7)) was added. 
The second central moment fi2 used for calculations was computed using first two {lf+mJ2(2)) and three 
{lf+mf2(3)) terms in its series expansion (j5.7ap . Finally, analytical estimates were completed by the 
third order {If+mf2(3)+mf3) VaR adjustment p.Sp . The estimation of the third central moment fi^ was 
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based on the first three terms of its series expansion listed in (|5.7bp . 

Comparison of the portfolio level results is presented in Table [U while VaR-based risk contributions 
on facility level are compared in Figure [T] The results presented in Figure [T] are sorted in increasing from 
left to right order using a scalar product of the principal vector Y and factor loadings vectors /3 as a 
parameter. 



If lf+mf2(2) lf+mf2(3) If+mf2(3)+mf3 

concentrated -5.2% -0.9% -0.8% -0.1% 

diversified -1.5% -0.1% -0.1% 0.0% 

Table 1: Relative differences between analytical approximation and Monte 
Carlo simulation on portfolio level. 

The following conclusions can be drawn based on the results of the numerical tests. Overall, the 
analytical approximation produces excellent results. On a portfolio level, a slight underestimation of 
VaR (economic capital) by a single factor approximation is observed for concentrated portfolios. The 
situation is improved by higher order corrections, whose contributions lead to very precise results. In case 
of VaR contributions, the higher order corrections to the single factor approximation (second and third 
order VaR adjustments) are are necessary to achieve high accuracy. The resulting analytical estimates 
of the VaR contributions are just 1-2% different from the Monte Carlo based estimates. 

5 Idiosyncratic risk 

The asymptotic multi-factor framework described in the previous sections was built on the conditional 
expectation series expansion. As a result, the idiosyncratic risk of the portfolio has been wiped out 
(averaged over) and portfolio risk measures were expressed in terms of the systematic components. In 
this section it is shown how the framework described so far can be extended to cover the idiosyncratic 
risk components. 

5.1 Idiosyncratic contributions 

Using the notations of the previous section and introducing the idiosyncratic value component Vga , the 
full portfolio value V can be written as 

where 

Vifim) = {V{vi,V*,0),.,^ (5-2a) 

Vmf{v*\m) = {v{m,v*,0)^-{vim,v*,0)^>.,i (5.2b) 

Vga{am,v*) = y{m,^\i)~{y{m,^\i))^■ (5.2c) 

The (. . .) in the above stands for average. Also, by construction, 

(V^<,a)5 = 0, (K„/),„. =0. (5.3) 



Following iPvkhtinl 's (|2004l) approach, one can treat the sum Vmj + Kja as a small correction to Vi/ 
and apply the results of Section |31 The VaR and ES contributions arc then expressed in terms of V\f 
and central moments of Vmf + ^ga- For the second and the third central moments considered here, one 
can write the following 

M2[K™/+K,a] = /^2[Ku/] + (M2[K,a(C)]),,. (5.4a) 

M3[K„/ + Vga] = /^3[K„/] + 3 (Kn/ . ^2 [K,a (0]),* + (M3 [^ga (C)]) ,. . (5.4b) 
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The first of tlie above is a well-known l ow of total varia nce, while the second one follows from a more 
general low of total cumulance (see, e.g.. iBrillingeil Il969l) . 



Thus, to compute the idiosyncratic contribution to the second and the third VaR and ES adjustments, 
one needs to compute the second and the third conditional central moments of the idiosyncratic compo- 
nent Vga- Only these moments, not a detailed information abut idiosyncratic components, are needed to 
complete the task. 

5.2 Conditional idiosyncratic moments 

Taking advantage of the conditional (on systematic factors) independence of the idiosyncratic components, 
the second and the third idiosyncratic central moments can be computed as a simple sum of individual 
(from each loan in the portfolio) contributions. The conditional expectation series expansion (|4.10[) can 
be applied not just to the value , but to its powers: 



W(r?) = 5]^z;rUft)fc,(A)fe.-.-He^^2...(^^) (5.5^) 



^K'y) = EtT^^ (/''Ofc.(Ak---He^'^2...(^^,) (5.5b) 



^Kv) = 5^^«r(ft)fe.(/3.)fe....He^^^2...(^^)^ (5.5e) 



where 



^f)= /«,(6)He„(e)^d6, u,^'' ^ j vKm^nief-^^^, ^^^ = /«f(6)He„(e)^d6 (5.6) 



The above expansions can be used to compute the conditional central moments of the ith loan in the 
portfolio 



((M2).),. = \«f;,,. -(^■^>r (5-7a) 



((M3).)„. = {vl)^^^2,{vl-v^)^^+2{vl-vl-v-)^^, (5.7b) 

and the mixed term in (j5.4bp 

3 {Vmf ■ (M2)«),. = 3 IVmf ■ ^)^^ ~i{VmfVl- W),. (5.8) 

The averages (. . .) , can be calculated using the techniques developed in Section 2] For example, 

an 

V^i.(«'He„(r/i) (5.9a) 

n 

E^! (E 0-^"^He„_.(,O(/3i)"-'= (V^^^)! 
f; Q«;(")He„_fe(r^i)(/?i)"-^ (\/l^) 



n! 
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fc,Z,m \ 2 )■ \ 2 )■ \ 2 )■ 



(Vi ■ Vi ■ Vi 



00 

' n 



^^^Y. L )«^"^Hc„_,(ryi)(/3i)"-'= ( y'l - PI 
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The central moments fi2 and /is of the portfoUo are calculated by summing over the ((/i2)i)„. and 
((/^3)j),i* i^ (|5.7a[) - (j5.7bp . The derivatives of these central moments are trivially calculated in a similar 
fashion using the property Hc„(77i) = n-He„_i(77i). Once the central moments and their derivatives have 
been calculated, the VaR and ES adjustments are obtained using p.7p - (|3.8p and (|3.10[) - p.lip . 

Due to the conditional independence, the idiosyncratic contribution is a sum of contributions from 
individual loans; hence, the amount of calculations needed is linear in a number of loans in the portfolio. 

5.3 Idiosyncratic risk allocation 

Using the techniques presented so far, allocation of idiosyncratic risk is a straightforward (although 
somewhat laborious) task. The Euler principle 12.31 can be applied to the idiosyncratic risk component in 
a way similar to the one described in Section 14.41 

Without going into the details of the calculations, let us emphasize the following peculiarity which 
never receive enough attention in the literature before. As a consequence of the conditional independence 
concept used, the portfolio level idiosyncratic risk contribution is calculated as a sum of the individual 
contributions from the underlying loans. These individual contributions, however, arc not equal to the risk 
contributions (referred to as marginal risk contributions) calculated using the Euler allocation principle. 
This is due to the fact that the partial derivative in Wi-^^ is applied not only to the central moments 
(and their derivatives), but also to the V',V",... terms in (j3.7|) - (|3.8p and p.lOp - p.lip . For example, 
a loan with pi = 1 in (|2.ip . i.e. with the value depending on systematic risk factors only, gives a zero 
contribution to the portfolio idiosyncratic moments ^2 and fi^. Yet, its contribution may be significant 
due to the effect mentioned above. 

Another interesting empirical property of the marginal idiosyncratic contributions is worth being 
highlighteqj. Although the overall contribution of the idiosyncratic risk component is positive and not 
significant for moderately concentrated portfolios, the marginal idiosyncratic contributions may be both 
positive and negative. It is common for most loans in the portfolio to receive small (compared to sys- 
tematic part) negative corrections, while a few big exposures receive relatively high positive ones. This 
counterintuitive behavior is demonstrated in the next section. 

5.4 Numerical results 

The analysis presented here mainly focuses on the effects due to the idiosyncratic risk and its concen- 
trations. In contrast with Section 14.51 realistic portfolio was used to benchmark the proposed analytical 
techniques against Monte Carlo simulations. As before, the portfolio VaR and the marginal VaR contri- 
butions were the risk measures of interest. 

The portfolio consisted of 2,000 loans to distinct customers randomly selected from a loan portfolio 
of a large European bank. The set of common systematic factors covering 45 geographic regions and 61 
regions, as well as the yaluation functio n at horizon Vi{e) used in the experiment were similar to those of 



the PortfolioManager (jKealhoferl . 120011) model. 

Both the portfolio VaR and VaR contributions were estimated using unbiased Monte Carlo simulations. 
The confidence level was set at 99.9%. To reduce simulation noise, the VaR contributions were estimated 
as value expectations in the interval 99.85%-99.95%. The simulations consisted of 10^° (ten billion) 
scenarios. 

The analytical estimates used for comparison with the simulation-based ones were constructed as 
follows. The systematic part of VaR (syst) was used as a starting point. These estimates covered up 
to the third order VaR adjustment as described in Section 14.51 The second and the third estimates 
were computed by consequently adding the second {ga2) and the third (gaS) ordei|j idiosyncratic VaR 
adjustments to the systematic part. 

The relative differences in the portfolio VaRs are presented in Table[2] One can observe, that the initial 
underestimation of VaR by the analytical approximation of systematic risk only is improved significantly 



^This is not a generic property of any portfolio. However, the behavior of the marginal contributions described here can 
be observed in most realistic portfolios. 

®The mixed term in Il5.4bt was included in the third order idiosyncratic contribution. This term gave negligible contri- 
bution in the test presented. It is, however, not clear if the term can be neglected in a general case. 
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Figure 2: Relative differences between Monte Carlo and analytical esti- 
mates of the VaR-based risk contributions. 
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Tabic 2; Relative differences between analytical and simulation-based es- 
timates of the portfolio VaR. 



by taking into account the second and the third order idiosyncratic contributions. Even without the third 
order corrections, the portfoho-level results are accurate enough for any practical purposes. 

The situation is different if we compare the marginal VaR contributions. Figure [5] shows relative 
differences in the VaR contributions. The triangles correspond to the ten loans with the highest VaR 
contributions. These are the loans with the highest single name risk concentration in the portfolio. The 
top scatter plot in Figure [2] shows that significant differences between the systematic (analytical) and the 
full (simulation-based) VaR contributions exist for the concentrated exposures. Both the second and the 
third idiosyncratic corrections need to be taken into account to obtain precise estimates. 

One can also observe the following property of the VaR contributions mentioned before. While a few 
biggest loans (in terms of VaR contributions) receive significant positive contributions to their VaRs due 
to idiosyncratic risk, most loans receive small negative contributions. 
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Overall, the results are consistent with those of Section 14.51 Applying the second and the third VaR 
adjustments to both systematic and idiosyncratic risk components leads to excellent results. On the 
portfolio level the results may be considered as exact. Computation of the marginal VaR contributions 
is a more challenging task; however, very high accuracy is achieved by taking into account higher order 
(the third order in this case) VaR adjustments. 

6 Summary 

The analytical framework for stru ctured credit p ortfolio models presented here is an extension and im- 
provement of the one developed bv lPvkhtinI ( 20041 ). Second and third order VaR and ES adjustments were 



considered. The default-only case considered by Pykhtin was extended to the case of arbitrary valuation 
function at horizon. The problem of quadratic (in portfolio size) complexity of Pykhtin's multi-factor ad- 
justment has been solved. High accuracy of the proposed technique was demonstrated by benchmarking 
with Monte Carlo simulations. The realized performance of the analytical approximation allows it to be 
considered as not just a supplement, but a substitute to the conventional simulation-based calculations. 
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